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Mott Transition for Strongly-Interacting ID Bosons in a Shallow Periodic Potential 


G. Boeris,* L. Gori,^ M. D. Hoogerland,^ A. Kumar,^ E. Lucioni,^ L. Tanzi,^ M. Inguscio,^’"*^ 

T. Giamarchi,^ C. D’Errico,^ G. Carleo,* G. Modugno,^ and L. Sanchez-Palencia* 

^Laboratoire Charles Fabry, Institut d'Optique, CNRS, Univ Paris Sud 11, 

2 avenue Augustin Fresnel, F-91127 Palaiseau cedex, France 
^LENS and Dipartimento di Fisica e Astronomia, Universitd di Firenze, and CNR-INO 50019 Sesto Fiorentino Italy 
^Department of Physics, University of Auckland Private Bag 92019, Auckland, New Zealand 

^INRIM, 10135 Torino, Italy 

^Department of Quantum Matter Physics, University of Geneva, 24 quai Ernest-Ansermet, 1211 Geneva, Switzerland 

(Dated: October 6, 2015) 

We investigate the superfluid-insulator transition of one-dimensional interacting Bosons in both deep and 
shallow periodic potentials. We compare a theoretical analysis based on Monte-Carlo simulations in continuum 
space and Luttinger liquid approach with experiments on ultracold atoms with tunable interactions and optical 
lattice depth. Experiments and theory are in excellent agreement. It provides a quantitative determination of the 
critical parameter for the Mott transition and defines the regime of validity of widely used approximate models, 
namely the Bose-Hubbard and sine-Gordon models. 


Introduction .— The interplay of repulsive interactions 
and a periodic potential in a quantum fluid triggers a 
superfluid-insulator transition known as the Mott transition, 
provided the potential period is commensurate to the inverse 
fluid density. The most familiar notion of Mott transition takes 
place in the limit of a deep periodic potential. In this case the 
lattice Hubbard model microscopically captures the dominant 
interaction and hopping processes, the strengths of which, U 
and J, strongly depend on the periodic potential amplitude V. 
Then, the Mott transition is driven by the competition of these 
sole two parameters at 7 ~ t/ [1, 2]. Quite strikingly in one 
dimension (ID) a Mott transition can exist even for a vanish¬ 
ingly small periodic potential provided the repulsive interac¬ 
tions are strong enough [3-6]. In the limit of a shallow poten¬ 
tial, its amplitude V becomes sub-relevant and the transition 
is mostly controlled by the interaction strength g alone [6]. 

Ultracold atoms provide a remarkable laboratory to study 
this physics [7, 8]. So far the Mott transition has been ob¬ 
served in both deep [9-13] and shallow [14] optical lattices. 
However, the characterization of the Mott transition in shal¬ 
low potentials remains a formidable challenge for both theory 
and experiments, with direct consequences not only in the ul¬ 
tracold atom realm but also in condensed matter for problems 
such as spin chains for instance [6, 15, 16]. 

On the theoretical side, while the Hubbard limit is now well 
documented [7, 8, 17] and its Mott transition has been exten¬ 
sively studied [18-24] its regime of validity beyond the deep- 
lattice limit is still unknown in ID, with ab-initio results hav¬ 
ing been reported so-far only for 3 dimensions (3D) [25]. In 
the limit of a vanishing potential, an estimate of the transition 
values may be found in the sine-Gordon model whose coeffi¬ 
cients are determined perturbatively [14, 26]. This, however, 
ignores the unavoidable renormalization of the field-theoretic 
coupling parameters by the potential, which may significantly 
affect the transition. On the experimental side, the Mott transi¬ 
tion has been clearly observed in the shallow lattice limit using 
modulation spectroscopy and transport measurements [14]. 
However, the experimental uncertainties did not allow to a 


precise determination of the phase diagram. 

In this Letter, we report the first quantitative joint theoret¬ 
ical and experimental investigation of the Mott transition for 
strongly-interacting ID Bosons in a shallow periodic poten¬ 
tial. Using continuous-space quantum Monte Carlo calcula¬ 
tions, we determine the exact quantum phase diagram. Our 
calculations confirm the field-theoretical universal predictions 
and provide, in addition, an accurate quantitative values of the 
critical parameters of the Mott transition. Experimentally, we 
perform transport measurements on a Bose gas with tunable 
interactions down to the limit of very shallow lattices and we 
analyze them with a phase slip based model to determine ac¬ 
curately the Mott transition. The numerical and experimental 
results are in excellent agreement and show significant devia¬ 
tion from the perturbative sine-Gordon theory using bare Lut¬ 
tinger parameters. 

Model and theoretical approach .— We consider zero- 
temperature interacting ID Bosons of mass m with a con¬ 
tact interaction of strength g, subjected to a periodic poten¬ 
tial y(x) =V sin^(kx) of spacing a = n/k and amplitude V. 
Both the large V and small V limits have the possibility of 
a Mott transition when the interactions are increased [17]. 
In spite of their qualitative different natures the two limiting 
cases are, however, expected to belong to the same univer¬ 
sality class for they both lead to the same low-energy sine- 
Gordon model [4, 6, 17, 27]. Within the Tomonaga-Luttinger 
liquid (TEL) approach, the homogeneous superfluid is param¬ 
eterized by the Luttinger parameter K, which characterizes 
the interaction strength. Eor weak interactions, the periodic 
potential is essentially irrelevant, except in renormalizing the 
effective value of the Luttinger parameter. Eor strong interac¬ 
tions, the TEL may be unstable upon introducing a periodic 
potential, which signals the Mott insulator phase. More pre¬ 
cisely, the Mott transition may be triggered either by changing 
the fluid density to commensurability at sufficiently strong in¬ 
teractions (Mott-5 transition) or by increasing the interactions 
at commensurability (Mott-U transition). The TEL theory 
predicts the universal critical values K], = 1/p^ and K,, = 2/p^ 
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Figure 1. (Color online) (a) QMC results for the number of particles per site, na, the compressibility K, the superfluid fraction /„ and the 
Luttinger parameter K, versus the chemical potential /i, for a periodic potential amplitude V = 2E^ and an interaction strength g = IfP' jma. 
The various curves correspond to increasing system sizes, Lla = 30, 50, and 100 (blue, green, and red respectively). The vertical dotted lines 
show the transition points determined from the criterion A), = 1. (b) Phase diagram in the g-fl plane for the fixed amplitude of the periodic 
potential V = 2E,. The black points joint by lines are determined from the = 1 criterion, while the colored points show the results found 
from the crossing point of the compressibility (red), the cusp of the compressibility (green), and the crossing point of the Luttinger parameter 
(blue). The red crosses are the hard-core limits. The inside of the lobe (red region) and the dashed line correspond a density of one particle per 
potential spacing. Inset: Mott gap versus interaction strength. The black dashed line is a fit to the BKT prediction (see text). The blue dashed 
line is the prediction of the BH model and the blue cross is the corresponding tip of the lobe. 


for the Mott-5 and Mott-t/ transitions respectively, where p is 
the commensurability order [5, 6 , 27]. 

However the TLL theory involves effective parameters that 
are not easily related to the microscopic Hamiltonian parame¬ 
ters and the critical curve is presently not quantitatively 
known. To precisely determine the Mott transitions, we use 
quantum Monte-Carlo (QMC) simulations. This allows us 
(i) to determine quantitatively the phase diagram in terms of 
the microscopic parameters and (ii) to compute explicitly the 
Luttinger parameter AT as a function of the microscopic ones 
and make the link with field theory. We use the same imple¬ 
mentation of the continuous-space worm algorithm [28, 29] 
in the grand-canonical ensemble as used in Ref. [30], which 
is numerically exact for all the physical quantities we study in 
the following [31]. 

Incommensurate transition .— We start with the incom¬ 
mensurate (Mott-5) transition, which may be triggered by 
changing the chemical potential ji. In order to accurately de¬ 
termine the critical point, several quantities are examined. The 
particle density n is computed directly in QMC and the com¬ 
pressibility K= dn/dji is computed independently from par¬ 
ticle number fluctuations. The hydrodynamic superfluid den¬ 
sity n, is found from the superfluid stiffness T„ defined as the 
response of the system to a twist of the boundary conditions, 
which is computed using the winding-number estimator [31]. 
We then deduce the superfluid fraction f^ = njn and the Lut¬ 
tinger parameter K = TZy/njc. 

The QMC results are shown versus the chemical poten¬ 
tial in Fig. 1 for y = 2£’, where £, = tP'k^/2m is the recoil 
energy, g = Ifi^/ma, and various system sizes. The density 
[Fig. l(al)] increases monotonically with the chemical poten¬ 


tial p. and exhibits a plateau at commensurability, na = \, 
where the superfluid density drops to zero. This is the signa¬ 
ture of the Mott-5 transition. The critical chemical potentials 
corresponding to the two edges of the plateau are accu¬ 
rately determined from the crossing points of the compress¬ 
ibility for different system sizes [Fig. I(a2)]. They can also be 
found from the drop of the superfluid fraction /„ which yields 
similar values for pf" [Fig. I(a3)]. At the Mott-5 transition, 
the Luttinger parameter is expected to exhibit the universal 
discontinuity from AT = 1 to AT = 0. Our data [Fig. I(a4)] are 
perfectly compatible with this prediction of the TLL theory. 

Repeating the same calculations for various values of the 
interaction strength, we find the Mott lobe in the g — p plane 
shown in Fig. 1(b). The black points and joining lines are de¬ 
termined from the AT^ = 1 criterion while the colored points are 
extracted from the crossing point of the compressibility (red), 
the cusp of the compressibility (green), and the crossing point 
of the Luttinger parameter (blue). The different methods yield 
results in excellent agreement all along the lobe within a few 
percents. The Mott gap A = — p^ is shown in the Inset of 

Fig. 1 versus the interaction strength. The black dashed line 
is a ht with the exponential closing A{g) oc exp(—h/y'g — 
predicted by the Berezinskii-Kosterlitz-Thouless (BKT) the¬ 
ory [32-34], with b and g^ as htting parameters. We hnd that 
the fit is good for all values of g except very close to the hard¬ 
core limit, g + 00 . The fit also yields a fair estimate of the 
critical value at commensurability, g^/(h^/ma) = 2.2 ± 0 . 1 , 
which corresponds to the tip of the lobe. 

Commensurate transition .— In order to characterize the 
commensurate (Mott-t/) transition more precisely, we now 
vary the interaction strength 7 = mg/H^n at commensura- 
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Figure 2. (Color online) Analysis of the correlation function gi {x) = 
() for VIE, = 2. The decay is exponential in the insu¬ 
lator (left inset, y =1) and algebraic in the superfluid (right inset, 
7 = 1.25). Left: Correlation length in the insulator. Right: Luttinger 
parameter in the superfluid extracted from the decay of the correla¬ 
tion function (black) and from the formula K = (red). The 

QMC results are the points, connected by straight lines to guide the 
eye. Note that the inflection of the correlation length curve close to 
the transition is due to finite size effects. 


Figure 3. (Color online) Phase diagram in the g-V plane at commen- 
surability (na = 1). The black and green points, both with errorbars) 
are the Monte Carlo and experimental results, respectively. The blue 
line is the BIT prediction using exact Wannier functions and the criti¬ 
cal value {J/U)c = 0.297 ± 0.01 (the shaded area corresponds to this 
error bar). The red cross indicates the critical interaction strength 
gc such that = 2 in the Lieb-Liniger model. Along the whole 
transition line the effective Luttinger parameter remains K = 2. The 
red line is the result of the bare sine-Gordon theory [26]. 


bility, na = 1. We compute the one-body correlation func¬ 
tion gi(x) = {\[f\x)\jf{0)), where \j/ is the field operator, 
from the statistics of endpoints on open world lines in the 
QMC [28, 29]. When increasing 7 along the line with na = 1 
(dashed line in Fig. 1), we observe a clear change of behavior 
of the gi function from algebraic for 7 < 7 , to exponential for 
7 > 7 (see Insets of Fig. 2). This is the signature of the Mott- 
U transition. The finite correlation length f in the insulating 
phase in shown in the left panel of Fig. 2. It is of only a few 
lattice sites long for strong interactions and increases up to a 
value comparable to the system size for 7 ~ 2. This is com¬ 
patible with the expected divergence of the correlation length 
at the transition. In the superfluid phase, the algebraic decay 
of the correlation function is compatible with the TLL theory 
prediction gi(x) 0 = 1/x*/^^. The two values of the Luttinger 
parameter found from a fit to this prediction and from the ther¬ 
modynamic prediction K — n^/rifK are in good agreement (see 
right panel of Fig. 2). When increasing the interaction towards 
the insulating phase, the Luttinger parameter decreases down 
to /T ~ 2 as predicted by the TLL theory. 

To locate the Mott-t/ transition point accurately, we resort 
to the BKT renormalization group equations to perform the 
finite size scaling of the Luttinger parameter [31]. The re¬ 
sults are shown in Fig. 3 (black points). In the strong potential 
limit, the results are compatible with the prediction of the BH 
model with the critical value {J/U)^ — 0.297 ±0.01 [21, 22] 
and the hopping J and interaction strength U calculated from 
the exact Wannier functions [20]. In the vanishing potential 
limit, the results converge to the critical value 7 (F = 0) ~ 3.5 
(red cross in Fig. 3) found from the exact relation K{g) for 
the integrable Lieb-Liniger model [35]. However, in the inter¬ 


mediate regime, we find a strong deviation from the pinning 
transition line (red dashed line) computed in Ref. [26] from 
the perturbative sine-Gordon theory using the bare Luttinger 
parameters. This shows that to quantitatively obtain the phase 
diagram the renormalization of the Luttinger parameters by 
even relatively weak interactions is significant and cannot be 
ignored in the perturbative field theory. 

Experiment .— Having quantitatively characterized the 
Mott transition from theory, we now turn to the experimental 
investigation. The experiment starts with a Bose-Einstein con¬ 
densate of 35 000 atoms with tunable scattering length at 
a broad Feshbach resonance [36]. The BEC is split into about 
1000 vertical ID tubes by adiabatically loading a strong hor¬ 
izontal 2D optical lattice. Each tube contains on average 36 
atoms and the transverse trapping frequency, Ox = 27r x 40 
kHz corresponds to an energy higher than all other energy 
scales, realizing an effective ID geometry. In the longitudi¬ 
nal direction we then adiabatically rise a weak optical lattice 
with spacing a = Xjl =532 nm and normalized amplitude 
V jE, ranging from 1.0(1) to 4.0(4). The temperature of the 
system, T ~ 30 nK, is below the ID degeneracy temperature 
To — 50nK [37]. A magnetic field along the vertical direction 
holds the system against gravity and a residual harmonic trap 
potential, with frequency (Oy = 2nx 160 Hz, makes it inho¬ 
mogeneous. By varying the 3D scattering length a^o, we can 
tune the Lieb-Liniger parameter 7 in the range 0.07 — 7.4. The 
system parameters are chosen to obtain a mean atom number 
per potential period {na) =1. This implies that in most of the 
tubes there are one or two regions with local commensurate 
density, na=l, which can undergo a Mott-t/ transition. There 
is however a fraction of the tubes for which na<\ which can- 
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Figure 4. (Color online) Critical momentum as a function of the 
scattering length for four periodic potential depth values: V/E,= l 
(red), 2 (orange), 2.8 (green), 4 (blue). A piecewise second-order 
polynomial fit determines the critical values for the SF-MI transition: 
respectively ac/flo=392± 12, 214±6, 160±30, 122±8. Inset: time 
evolution of the momentum distribution peak p for =109ao and 
V /E^=2. The solid line is the theoretical damped oscillation fitting 
the data p < p^ before the dynamical instability sets in. The error 
bars represent the root mean square of the imaging resolution and 
statistical uncertainties. 


not become insulating [31]. 

To detect the Mott transition we excite a sloshing motion 
of the system through a shift of the trapping potential, ob¬ 
tained by suddenly switching off the magnetic field gradi¬ 
ent [14, 38]. We let the atoms evolve in the trap for a vari¬ 
able time f, after which all optical potentials are switched off 
and time-of-flight absorption images are recorded. An ex¬ 
ample of the time evolution of the momentum distribution 
peak p{t) is shown in the inset of Fig. 4: One can notice 
an initial increase of p followed by a subsequent decrease. 
We analyze this behavior in the frame of a phase slip based 
model [38, 39]. Phase slips, i.e. the dominant excitations in 
ID, make the system dynamics dissipative already at small 
momenta: p{t) can be fit with a damped oscillation function 
of the form p{t) = sin(a)'f), where the frequency (o' 

is renormalized by the damping rate G and by the lattice [31]. 
At larger momenta one observes a clear deviation from such 
behavior, with a sudden increase of the damping that signif¬ 
icantly reduces the growth of p. We identify the momentum 
where the experimental data points deviate with respect to the 
theoretical curve as the critical momentum p^ for the occur¬ 
rence of a dynamical instability, driven by a divergence of the 
phase slip rate [38]. The critical momentum p^ is expected to 
vanish at the superfluid-insulator transition [39]. 

The behavior of p^ as a function of the scattering length a,o 
is reported on Fig. 4 for several values of the lattice depth. The 
measured p^ initially decreases for increasing a^o, and then 
reaches a finite constant value. We interpret the onset of the 
plateau as the Mott transition for the commensurate regions 
of the system, which drives all the corresponding tubes into 


an effective insulating regime (i.e. transport along individ¬ 
ual tubes is globally suppressed). The fraction of tubes that 
does not reach the critical density na=l keeps instead mov¬ 
ing also beyond this point, originating the observed plateau 
for p^. This interpretation is corroborated by the observed in¬ 
crease of Pc at the plateau for decreasing V, since the increase 
of the interaction strength that is necessary to reach the insu¬ 
lating regime produces an overall decrease of the density of 
the ID systems, hence an increase of the fraction of tubes that 
does not reach na=l. For each set of measurements with a 
given value of V, we therefore identify the critical scattering 
length value to enter the insulator regime by determining 
the beginning of the plateau with a piecewise fit. We use a 
second-order polynomial fit, which is justified by the phase 
slip based model [31, 39, 40]. We clearly see that as V de¬ 
creases, Uc - and thus also jc - increases. 

For each periodic potential depth we get the Mott-t/ transi¬ 
tion point converting into 7 , for na = 1. The experimental 
results are shown as green points in Fig. 3. Within our uncer¬ 
tainties [41], the experiment is in very good agreement with 
the numerical simulations, except for the case VlE=\, where 
the finite size of the system might start to play a role. Our re¬ 
sults are also consistent with those reported in Ref. [14] within 
their uncertainties. Our experiment confirms the clear devia¬ 
tion of the transition line from the bare sine-Gordon prediction 
observed in the QMC results. Note that, surprisingly enough, 
while the BH model is justified only for V ^ E„ both numer¬ 
ics and experiments show that the BH prediction for the Mott- 
U transition is quite accurate down to the limit F —> 0. This 
agreement is however rather accidental and the breakdown of 
the BH model is manifest in other quantities. For instance, the 
BH prediction for the Mott gap deviates significantly from the 
exact QMC result [see Inset of Fig. 1(b)]. 

Conclusions .— We have studied, both theoretically and 
experimentally, the Mott transition of strongly-interacting ID 
Bosons in a periodic potential from deep to shallow potentials. 
Our ab-initio QMC calculations validate the field-theoretic 
universal predictions and in addition provide a quantitative de¬ 
termination of the phase diagram. It shows that the renormal¬ 
ization of the Luttinger parameter is significant even for weak 
periodic potentials. The numerical analysis give excellent 
agreement with experiments for the Mott-t/ transition. The 
experimental observation of our numerical phase diagram for 
the Mott-5 transition is still beyond reach for ultracold atomic 
systems due to the requirement of a fine control of atom num¬ 
ber in box-shaped potentials. In spite of recent progress in 
that direction [42, 43], it remains a great challenge for future 
studies. 

This research was supported by the EU FET-Proactive 
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(EP7/2007-2013 grant No. 256294), Marie Curie lEE 
(PP7/2007-2013 grant No. 327143), MIUR (grant 
No. RBPR12NLNA), and the Swiss NSE under Division II. 
Numerical calculations were performed using HPC resources 
from GENCI-CCRT/CINES (grant No. c2015056853) and 
the GMPCS cluster of the LUMAT federation (ER LUMAT 
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2764), and make use of the ALPS scheduler library and 
statistical analysis tools [44—46]. 

Note added .— During the completion of this manuscript, 
a preprint appeared reporting the numerical study of the Mott- 
U transition with results consistent with ours [47]. 
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-Supplemental Material- 

Mott Transition for Strongly-Interacting ID Bosons in a Shallow Periodic Potential 

In this Supplemental Material we give additional details on the methods and the analysis we have performed in the main 
Letter. In particular we provide a detailed description of i) the quantum Monte Carlo calculations, ii) the determination of the 
superfluid density, iii) the finite-size scaling analysis, and iv) the experimental methods. 


Quantum Monte Carlo calculations 


Superlluid density 


We consider ID interacting Bosons at zero temperature 
governed by the continous-space Hamiltonian, 


H = dx 


^ tP' e 

2m 2 


, (SI) 


where (//(jt) is the Bose held, m is the particle mass, g is the 
contact interaction strength, and V{x) = VmP'ikx) is a peri¬ 
odic potential of spacing a — n/k and amplitude V. 

We study the interacting quantum system from first prin¬ 
ciples, by means of quantum Monte Carlo (QMC) numerical 
simulations. Path-integral QMC approaches allow us to com¬ 
pute the exact thermodynamic properties of a generic, inter¬ 
acting bosonic quantum system in any dimensions. This is 
achieved by means of the path-integral representation of the 
grand-canonical partition function Z = where 

P = I/ksT with kg the Boltzman constant and T the temper¬ 
ature, ji is the chemical potential, and iV = /ck y/ is the 
total particle number operator. The quantum partition func¬ 
tion is expressed as an equivalent classical partition function 
of interacting polymers living in dimension D+\, where the 
additional dimension is the imaginary-time direction [1]. The 
equivalent classical partition function can be treated stochas¬ 
tically by means of Monte Carlo sampling. An efficient way 
of sampling the associated partition function is given by the 
worm algorithm [2, 3]. 

In our implementation, the imaginary-time propagator en¬ 
tering the path-integral representation is written as 


/=! 1=1 

r( 2 ) 




‘<J 


where = \Jm/2KTpe x e is the free- 

particle propagator. The one-body contribution V/ip) = 
/p£*^'(ri,r') is computed exactly by 
matrix squaring [1]. The two-body interaction is taken 
into account at the pair-product level, thanks to the explicit 


_^(2) /j. I* . I*' 

expression for the two-body propagator e ^ i’ = 

which is known [4]. The systematic error coming from the 
discretization of the path-integral along the imaginary-time 
direction is smaller than the statistical errorbars reported in 
our results. 


The stiffness Y = — 


niL d^XnZ 


measures the response of 


the system when the periodic boundary conditions are twisted 
by an angle Oq. It is directly computed by QMC using the 
usual winding number estimator [1]. However, this quantity 
strongly depends on the ratio j3/L [5]. The relevant quantity 
to study is the superfluid density n, appearing in the hydrody¬ 
namic action 


S[0] = 


rL f-np 

/ dx 

/ dT 

JO 

'0 


fPn^ 




2m 2 


(S3) 


Computing the partition function Z = J cor¬ 

rectly taking into account winding configurations of the form 
9{x, t) = 2%kxlL with k gZ, one deduces the relation 


T, = n J 1 - 4-n^n, 


mL g—lTt^n^k^Tp^/mL 


(S4) 


^k=-o 


In our calculations, we invert this relation and extract n, from 
the stiffness computed by QMC. 


Finite size scaling 


To precisely determine the Mott-1/ critical interaction 
strength we study the evolution of the Luttinger parame¬ 
ter K for the increasing system sizes Lla = 30, 50, and 100. 
When L is increased, the temperature T is lowered by keep¬ 
ing the ratio j3 jL constant in order to consistently compute 
ground-state properties. We resort to the known Berezinkii- 
Kosterlitz-Thouless (BKT) renormalization group (RG) equa¬ 
tions to perform the finite size scaling and determine the crit¬ 
ical interaction strength in the thermodynamic limit. In par¬ 
ticular, the finite-size behaviour of the Luttinger parameter K 
and of renormalized potential strength Y are governed by the 
renormalization equations [6-10], 


df 16 \EJ 


— = {2-K)V. (S5) 


Dividing these two equations one by another, one finds that 
the quantity 





y 

16 \ E, 


(S6) 


is conserved along the RG trajectories. The separatrix be¬ 
tween the superfluid and the insulating phases then corre¬ 
sponds Xo K = 2 and V = 0, i.e. C = 1- In practice, we use 
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the values the Luttinger parameter K\ and K 2 for two system 
sizes Li and L 2 , and compute ^ using the integral relation 

The Mott-t/ critical point is then determined using the crite¬ 
rion C=l. 


Experimental methods 

In the experiment, a Bose-Einstein condensate of atoms 
is split in about 10^ ID subsystems (potential tubes) by means 
of a 2D optical lattice. The overall Thomas-Fermi distribu¬ 
tion of atoms in the tube labelled by i and j along the two 
horizontal directions is given by Nij = A^o.o[l ~ 2.7iNofi{i^ + 
where = 3.5 x lO'* is the total atom number 
and A^op — 50 is the atom number in the central tube. On A^r 
there is a 6 % statistical error while a 30% systematic uncer¬ 
tainty is due to the imaging calibration. Along the vertical 
direction of the tube, where we load the weak periodic po¬ 
tential, the trapping frequency is w, = 2n x 160 Hz. A good 
estimate of the mean atomic density for each tube is provided 
by the largest of the Thomas-Fermi and the Tonks value [11]. 
The mean site occupation {no) is then calculated by averag¬ 
ing over all the tubes. During the lattice loading we employ an 
optimal value of the 3D scattering length to get commensura- 
bility on average, i.e. {no) = 1. Given (O^ and At, the optimal 
value is 030 = 220 ao- 

After the lattice loading, by varying the 3D scattering 
length fl 3 D, we tune the ID scattering length aio = £^(1 — 
1.03a3D/f±)/2a3D, and thus the Lieb-Liniger parameter 7 , 
where = \/h/mC 02 _ is the harmonic transverse size of a 
tube, fixed by the 2D lattice. The value of 7 is estimated 
by averaging over all the tubes 1 /{nmaxaciw), where the peak 
occupation n„axCi is the largest of the Thomas-Fermi and the 
Tonks value. The mean peak occupation, i.e. the peak occu¬ 
pation averaged over all the tubes, is ~ 1.2 and the peak 
density in the central tube is less than two, thus preventing 
localization mechanisms at occupations other than one. De¬ 
spite the inhomogeneity of our system, at sufficiently strong 
interactions we clearly observe a suppression of the system 
dynamics. We interpret the latter as due to the fact that within 
each tube, a part of the atoms reaches the localization con¬ 
dition na = l stopping also the remaining adjacent parts with 
different occupation. We estimate that about one quarter of the 
atoms resides in tubes where the occupation is always na <\, 
justifying the plateau in shown in Fig. 4 of the main paper. 

The complete model we use to fit the time dependence of 
the momentum, pit), as shown in Fig. 4 of the main paper, is 
p{t) = p„„e^^^sin(a)'f) with amplitude = m*o}*'^zo/0)', 
frequency (o’ = Vand damping rate G. Here m* is 
the effective mass due to the shallow lattice, zo — 3 pm is the 
trap displacement, and co* = CO^^/m jm* is the lattice renor¬ 
malized frequency. 


In the absence of an exact theoretical model for the criti¬ 
cal momentum p^ to enter the dynamical instability regime at 
finite interaction, we use a quantum phase slips based model 
to predict the interaction dependence of p^. Along the lines 
of Ref. [12] we use the equation for the quantum phase slip 
nucleation rate [13, 14] 


r =BL{JJ)\/naJU \/ k/ 2 — pX/2hx 


\l.\{Kl2-pXl2h)^l^ 


Y 2n^lJ jnaJ 

exp —7.1 \/ naJ/U {7z/2 — pX /2/i)^/^ 


(S 8 ) 


Here L{U) ~ 2t/*/^ is the mean length of the tubes and B is 
a phenomenological constant. Assuming that the system en¬ 
ters the strongly dissipative regime when the nucleation rate T 
exceeds a constant critical value, from Eq. (S 8 ) we can obtain 
p^ for each given value of the interaction U. Since this re¬ 
lation is valid in the Bose-Hubbard regime, we can employ 
it to estimate the critical momentum Pc{U) for the experi¬ 
mental case with V/E^ = 4. The phenomenological param¬ 
eter B is arbitrarly adjusted to reproduce the measured p^ at 
one interaction value, U = 2.47. In Fig. SI we compare the 
measured and the predicted p^ as a function of the scattering 
length a 3 D U. The quantum phase slip-based model well 
reproduces the experimental behaviour of p^ with the interac¬ 
tion, showing a quadratic dependence of on 030 . This result 
justifies the choice of a quadratic polynomial fit in Fig. 4 of 
the main paper, at least for the measurements with larger V. 
For shallower lattices the Tomonaga-Luttinger Liquid model 
for the nucleation rate could in principle be used [13], but this 
would require the Luttinger parameter A as a function of the 
interaction 7 in the presence of the lattice, which is still un¬ 
known. 
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Figure SI. For V/£r = 4 the experimental critical momentum (blue 
dots) is compared to the one predicted by the quantum phase slip 
based model, which applies up to the critical scattering length Oc to 
enter the MI regime (gray dashed line). For a > a^, Pc has instead a 
plateau due to tubes with occupation na < I (black dotted line). 




